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Abstract 

We consider wireless mesh networks and the problem of scheduling 
the links of a given set of routes under the assumption of a heavy-trafSc 
pattern. We assume some TDMA protocol provides a background of syn- 
chronized time slots and seek to schedule the routes' links to maximize 
the number of packets that get delivered to their destinations per time 
slot. Our approach is to construct an undirected graph G and to heuris- 
tically obtain node multicolorings for G that can be turned into efficient 
link schedules. In G each node represents a link to be scheduled and the 
edges are set up to represent every possible interference for any given set 
of interference assumptions. We present two multicoloring-based heuris- 
tics and study their performance through extensive simulations. One of 
the two heuristics is based on relaxing the notion of a node multicoloring 
by dynamically exploiting the availability of communication opportunities 
that would otherwise be wasted. We have found that, as a consequence, 
its performance is significantly superior to the other's. 

Keywords: Wireless mesh networks, Link scheduling, Node multicolor- 
ings. Scheduling by edge reversal. 
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1 Introduction 



Owing to their numerous advantages, wireless mesh networks (WMNs) consti- 
tute a promising solution for community networks and for providing last-mile 
connections to Internet users pQ [31 [H]- However, like all wireless networks 
WMNs suffer from the problem of decreased capacity as they become denser, 
since in this case attempting simultaneous transmissions causes interference to 
increase significantly |15l I25| . One common solution to reduce interference is 
to adopt some contention-free TDMA protocol [11] and schedule simultaneous 
transmissions for activation only if they do not interfere with one another. Do- 
ing this while maximizing some measure of network usage and guaranteeing 
that all links are given a fair treatment normally translates into a complicated 
optimization problem, one that unfortunately is NP-hard P|. 

This scheduling problem has been formulated in a great variety of manners 
and has received considerable attention in the literature. Prominent studies 
include some that seek to calculate the capacity of the network [HI [13] , others 
whose goal is the study of the time complexity associated with the resulting 
schedules [H], and still others that aim at scheduling transmissions in order to 
achieve as much of the network's capacity as possible [T0l[2l [28l[T2[[l6l[27l[23l[30l 
[29] . One common thread through most the latter is that, having adopted a graph 
representation of the network and of how the various transmissions can interfere 
with one another, a solution is sought through some form of graph coloring. 
More often than not the transmissions to be scheduled are represented by the 
graph's nodes and then node coloring, through the abstraction of an independent 
set to represent the transmissions that can take place simultaneously, is used. 
But sometimes it is the graph's edges that stand for transmissions, in which 
case edge coloring is used, building on the abstraction of matchings to represent 
simultaneity [3]. 

Here we consider a variation of the problem which, to the best of our knowl- 
edge, is novel both in its formulation and in the solution type we propose. We 
start by assuming a WMN comprising single-channel, single-radio nodes and 
for which a set of origin-to-destination routes has already been determined, and 
consider the following question. Should there be an infinite supply of packets at 
each origin to be delivered to the corresponding destination in the FIFO order, 
and should all nodes in the network be endowed with only a finite number of 
buff'ers for the temporary storage of in-transit packets, how can transmissions 
be scheduled to maximize the number of packets that get delivered to the desti- 
nations per TDMA slot without ever stalling a transmission, by lack of buffering 
space, whenever it is scheduled? This question addresses issues that lie at the 
core of successfully designing WMNs and their routing protocols, since it seeks 
to tackle the problem of transmission interference when the network is maxi- 
mally strained. The solution we propose is, like in so many of the approaches 
mentioned above, based on coloring a graph's nodes. Unlike them, however, we 
use node multicolorings instead p6^ , which are more general and for this reason 
allow for a more suitable formulation of the optimization problem to be solved. 
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Given the origin-to-destination routes (or paths) to be used, we begin in 
Section [2] with a precise definition of a schedule and a precise formulation of the 
problem. We also show, through an example, that had the problem been for- 
mulated for network-capacity maximization, a conflict with the requirement of 
finite buffering might arise. Then we move, in Section [31 to the specification of 
the undirected graph that underlies our algorithm's operation. One assumption 
in that section, and therefore throughout most of the paper, is that the commu- 
nication and interference radii are the same for the WMN at hand. Moreover, 
we also assume that the tenets of the protocol-based interference model [51 [53] , 
including the possibility of bidirectional communication in each transmission, 
are in effect. In Section [?] we guide the reader through various multicolor- 
ing possibilities, which culminates in Section [5] with a preliminary method for 
scheduling, borrowed from the field of resource sharing [5] . Improving on this 
preliminary method with the goals of the problem formulated in Section [2] in 
mind finally yields our proposal in Section [51 This proposal, essentially, stems 
from a slight relaxation of the notion of a node multicoloring. The subsequent 
two sections are dedicated to the presentation of computational results, with 
the methodology laid down in Section [71 and the results proper in Section [51 
Discussion follows in Section [51 and we close in Section [TUl 

2 Problem formulation 

We consider a collection Vi,V2, ■ ■ ■ jPp of simple directed paths (i.e., directed 
paths that visit no node twice), each having at least two nodes (a source and 
a destination). These paths' sets of nodes are Xi, X2, . . . , Xp, respectively, not 
necessarily disjoint from one another, and we let X = Up^^ Xp. Their sets of 
edges are Yi,Y2, . . . ,Yp and we assume that, for p ^ g, a member of Yp and one 
of Yq are distinguishable from each other even if they join the same two nodes 
in the same direction. Letting Y = lj^=i then see that Y may contain 

more than one edge joining the same two nodes in the same direction (parallel 
edges) or in opposing directions (antiparallel edges). 

Our discussion begins with the definition of the directed multigraph D = 
{X, Y) , where all P directed paths are represented without sharing any directed 
edges among them. An example is shown in Figure [U We take D to be repre- 
sentative of a wireless network operating under some TDMA protocol. In this 
network, each of paths 1^1,1^2, ■■■ , Vp is to transmit an unbounded sequence of 
packets from its source to its destination. Such transmissions are to occur with- 
out contention, meaning that whenever an edge is scheduled to transmit in a 
given time slot no other edge that can possibly interfere with that transmission 
is to be scheduled at the same time slot. We assume that each transmission 
sends at most one packet across the edge in question (more specifically, it sends 
exactly one packet if there is at least one to be sent but does nothing otherwise) . 
We also assume that each transmission may involve the need for bidirectional 
communication for error control. 
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Figure 1: A set of P = 3 directed paths (a) and the resulting directed multigraph 
D{h). 



We call a schedule any finite sequence S = {Sq, Si, . . . , Sl-i) such that 
Se C Y for < £ < L - 1, provided U^J„ Si = Y and moreover no two 
concurrent transmissions on edges of the same Si can interfere with each other. 
To schedule the transmissions according to S is to cycle through the edge sets 
^o, Si, ... , Sl-i, indefinitely and in this order, letting all edges in the same set 
transmit in the same time slot whenever that set is reached along the cycling. 
Given S, we let length(<S) = L and denote by delivered(iS) the number of packets 
that can get delivered to all paths' destinations during a single repetition of S 
in the long run (i.e., in the limit as the number of repetitions grows without 
bound). Of course, delivered(<S) is bounded from above by the number of times 
the P paths' terminal edges (those leading directly to a destination node) appear 
in S altogether. 

Before we use these two quantities to define the optimization problem of 

finding a suitable schedule for D, we must recognize that our focus on the 
source-to-destination packet flows on the paths Vi,V2, - ■ ■ ^Vp carries with it 
the inherent constraint that the nodes' capacity to buffer in-transit packets 
cannot be allowed to grow unbounded. We then adopt an upper bound B on the 
number of in-transit packets that a node can store for each of the paths (at most 
P) that go through it. However, there is still a decision to be made regarding 
the effect of such a bound on the transmission of packets. One possibility would 
be to impose that, when it is an edge's turn to transmit it does so if and only if 
there is a packet to transmit and, moreover, there is room to store that packet 
if it is received as an in-transit packet. Another possibility, one that seeks to 
never stall a transmission by lack of a buffer to store the packet at the next 
intermediate node, is to only admit schedules that automatically rule out the 
occurrence of such a transmission. We adopt the latter alternative. 

The following, then, is how we formulate our scheduling problem on D. Find 
a schedule <S that maximizes the throughput 

^^^> - length(5) ' 

subject to the following two constraints: 

CI. Every node can store up to B in-transit packets for each of the source-to- 
destination paths that go through it. 
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Figure 2: A set of P = 3 directed paths (a) and the resulting directed muhigraph 
D (b). Using the schedule S such that Sq = {a,/}, 5*1 = {c, d}, 5*2 — {b}, 
S3 = {e}, 5*4 = {a,g}, and S5 = {h} causes unbounded packet accumulation at 
node 2 when constraint C2 is in effect, thus violating constraint CI. Enforcing 
constraint CI for some value of B causes constraint C2 to be violated. 



C2. Whenever an edge is scheduled for transmission in a time slot and a packet 
is available to be transmitted, if the edge is not the last one on its source- 
to-destination path then there has to be room for the packet to be stored 
after it is transmitted. 



2.1 Scheduling for maximum network usage 

Before proceeding, recall that, as mentioned in Section [1] the most commonly 
solved problem regarding the selection of a schedule S is not the one we just 
posed, but rather the problem of maximizing network usage. In terms of our 
notation, this problem requires that we find a schedule that maximizes 

without any constraints other than those that already participate in the defini- 
tion of a schedule. 

It is a simple matter to verify that solutions to this problem often fail to re- 
spect constraints CI and C2 of our formulation. This is exemplified in FigureO 



3 Graph transformation 

We wish to address the problem of optimizing T{S) exclusively in terms of 
some underlying graph. Clearly, though, the directed multigraph D is not a 
good candidate for this, since it does not embody any representation of how 
concurrent transmissions on its edges can interfere with one another. Our first 
step is then to transform D into some more suitable entity, which will be the 
undirected graph G — {N, E) defined as follows: 

1. The node set of G is the edge set Y of D. In other words, G has a 
node for every edge of D. Since Z) is a multigraph, a same pair of nodes 
«, j € X such that (i, j) G F or (j, i) e F may appear more than once as 
a member of N . 
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2. The edge set of G is obtained along the following four steps: 

i. Enlarge N by including in it all node pairs of D that do not cor- 
respond to edges on any of the P source-to-destination paths but 
nevertheless reflect that each node in the pair is within the interfer- 
ence radius of the other. We refer to these extra members of N as 
temporary nodes. 

ii. Connect any two nodes in N by an edge if, when regarded as node 
pairs from D, they share at least one of the nodes oi D. In other 
words, if each of the two pairs i,j £ X and k,l d X corresponds to 
a node of G (by virtue of either constituting an edge of D or being a 
temporary node), then the two get connected by an edge in G if at 
least one oi i — k, i — I, j — k, or j — I holds. 

iii. Connect any two nodes in N by an edge if, after the previous steps, 
the distance between them is 2. 

iv. Eliminate all temporary nodes from N and all edges from E that 
touch them. 

Together, these four steps amount to using G to represent every possible 
interference that may arise under the assumptions of the protocol-based 
model when communication is bidirectional. Graph G is also known as 
a distance-2 graph relative to I? [1]. The entire transformation process, 
from the set of P paths through graph G, is illustrated in Figure |31 

It follows from this definition of G that any group of nodes corresponding 
to parallel or antiparallel edges in D are a clique (a completely connected sub- 
graph) of G. Similarly, every group of three consecutive edges on any of the 
paths 7^1, 7^2, ■ • ■ , Vp corresponds to a three- node clique in G. As we discuss in 
Section |9l these and other cliques are related to how large T{S) can be under 
one of the scheduling methods we introduce. 

It is also worth noting that Steps 1 and 2 above are easily adaptable to mod- 
ifications in any of the assumptions we made. These include the assumptions 
that the communication and interference radii are the same and that communi- 
cation is bidirectional. Changing assumptions would simply require us to adapt 
Steps 2.i through 2. iii accordingly. 

4 Multicoloring-based schedules 

Graph G allows us to rephrase the definition of a schedule as follows. We call 
a schedule any finite sequence S = {So, Si, . . . , Sl-i) such that Si C N for 
0<e<L-l, provided IJfco = N and moreover every Si is an independent 
set of G. The appearance of the notion of an independent set in this definition 
leads the way to a special class of schedules, namely those that can be identified 
with graph multicolorings [26j . 

For q > 1, a. q-coloring of the nodes of G is a mapping from N, the graph's 
set of nodes, to N'^, where N is the set of natural numbers, such that no two of a 
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Figure 3: The graph-transformation process. We start with the directed muhi- 
graph D (a), to which the node pair labeled x is added as a dashed line to 
indicate the existence of interference that is not internal to any of the initial P 
paths. Panel (b) contains the undirected graph G as it stands after Step 2.ii. 
Panels (c) and (d) show G past Steps 2.iii and 2.iv, respectively. In these two 
panels, dashed lines are used to represent the edges added in Step 2.iii. 

node's q colors are the same and besides none of them coincides with any one of 
any neighbor's q colors. Of course, the set of nodes receiving one particular color 
is an independent set. If p is the total number of colors needed to provide G 
with a g-coloring, then A'^ is covered by the p independent sets that correspond 
to colors and every node is a member of exactly q of these sets. Therefore, 
letting L = p and identifying each with the set of nodes receiving color £ 
implies that to every g'-coloring of the nodes of G there corresponds a schedule 
S. 

These multicoloring-derived schedules constitute a special case in the sense 
that every node of G can be found in exactly the same number of sets (q) out 
of the L sets that make up the schedule. Clearly, though, there are schedules 
that do not correspond to multicolorings. For now we concentrate on those that 
do and note that delivcrcd(iS) < Pq always holds (recall that P stands for the 
number of origin-to-destination paths). That is, the greatest number of packets 
that the P terminal edges of D can deliver during the L time slots of schedule 
<S is g' per terminal edge. These schedules can be further specialized, as follows. 
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4.1 Standard coloring 



When q = I every node of G receives exactly one color and lengtli(iS) = L > 
x{G), where x(G') is the least number of colors with which it is possible to 
provide G with a 1-coloring, known as the chromatic number of G. Using T^[S) 
to denote T{S) in this case, we have 

T\S) < j^, (3) 



4.2 Standard multicoloring 

Coloring G's nodes optimally in the previous case is minimizing the overall 
number of colors. This stems not only from the fact that q = 1, but more 
generally from the fact that q is fixed. We can then generalize and define x'^(G) 
to be the least number of colors with which it is possible to provide G with 
a g-coloring. Evidently, x{G) = X^{G) < X^{G) < • • •, so the question of 
multicoloring G"s nodes optimally when q is not fixed can no longer be viewed 
as that of minimizing the overall number of colors needed (as this would readily 
lead to g = 1 and x{G) colors). Instead, we look at how efficiently the overall 
number of colors is used, i.e., at what the value of q has to be so that x'^{G)/q 
is minimized. This gives rise to the multichromatic number of G, denoted by 
X*(G) and given by X*i.C!) = iniq>ix'^{G)/q. Because this infimum can be 
shown to be always attained, we use minimum instead and let q* be the value 
of q for which x*(G) = x'^' {G)/<f- 

Using a q-coloring for scheduling amounts to having length(iS) ~ L > x^(G). 
In this case, letting T*{S) stand for T{S) yields 

T*(S) < < ^ (A) 

^ X^iG) - x«*(G) x*{G)- ^ ' 



4.3 Interleaved multicoloring 

A special class of g-colorings is what we call interleaved g-colorings [71 El [31] . If 
i and j are two neighboring nodes of G, let < C2 < ■ • ■ < be the q colors 
assigned to node i by some g-coloring, and likewise let < < • • • < be 
those of node j. We say that this g-coloring is interleaved if and only if either 
c'l < cj < < < • • • < < or cj < c'l < < c| < • • • < < for 
all neighbors i and j. If we restrict ourselves to interleaved g-colorings, then 
similarly to what we did above we use Xhit(^) denote the least number of 
colors with which it is possible to provide G with an interleaved g-coloring, and 
similarly define the interleaved multichromatic number of G, denoted by x*nt (G) , 
to be Xhitl'-') ^ hifq>i Xi',it(^)/'?- Once again it is always possible to attain the 
infimum, so we may take q* to be the value of g for which Xint(^) = XintiG) / q* ■ 
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As for scheduling based on an interleaved g-coloring, it corresponds to having 
length(5) = L> xfnti'^)- before, we use T*^^,{S) in lieu of T{S) and obtain 

T* (S) < < - P (5) 



4.4 Discussion 

It is a well-known fact that 



X(G) - xUG) - X*{G)- ' ' 

The first inequality follows from the definition of xtntiG), considering that ev- 
ery 1-coloring is (trivially) interleaved. As for the second inequality, it fol- 
lows directly from the definition of X*{G)- By these inequalities, should all of 
Eqs. ©-([I]) hold with equalities, we would have 

T\S) < TUS) < T*{S). (7) 

Obtaining equalities in Eqs. ©-([S]), however, requires both that delivered(5) = 
Pq for g = 1 or q — q* , as the case may be, and that length(5) = x'^{G) with 
the same possibilities for q or length(5) = xfntiG)- 

While the combined requirements involve the exact solution of NP-hard 
problems (finding any of xiG), xtntiG), and X*{G) is NP-hard; cf., respec- 
tively, [17], [7], and [14^), the former one alone is always a property of schedules 
based on multicolorings when buffering availability is unbounded. To see that 
this is so, first recall that the definition of delivered(5) refers to a repetition of 
the whole schedule as far down in time as needed for any transient effects to 
have waned. So, given any of the P source-to-destination paths, we can prove 
that delivered(iS) — Pq by arguing inductively about what happens on such a 
path during that future repetition of S. The basis case in this induction is the 
first directed edge on the path (the one leading out of the source) . The prop- 
erty that every appearance of this edge does indeed transmit a packet follows 
trivially from the fact that the source has an endless supply of new packets to 
provide whenever needed. Assuming that this also happens to the next-to-last 
edge on the path (this is our induction hypothesis) immediately leads to the 
same conclusion regarding the last edge, the one on which delivered(5) is de- 
fined. To see this, let e be the last edge and e~ the next-to-last one. Because S 
is repeated indefinitely, every time slot t sufficiently far down in time in which 
e appears is the closing time slot of a window in which both e and e~ appear 
exactly q times each. By the induction hypothesis, it follows that at least one 
packet is guaranteed to exist for transmission through e at time slot t. 

Buffering availability, however, is not unbounded, so we must argue for its 
finiteness. We do this by recognizing another important property of multicoloring- 
based schedules, one that is related to constraints CI and C2 introduced earlier. 
Because every edge of D (node of G) appears the exact same number q of times 



9 




Figure 4: A set of P = 4 paths (a), with dashed hnes indicating all node pairs 
representing ofF-path interference. The resulting graph G is shown in panel 
(b). Depending on the schedule S it is possible to obtain equalities in all of 
Eqs. (jni)-®. The schedules that achieve this while implying strict inequali- 
ties in Eq. © are: S = {{a,d,g},{bj,h},{c,e}) for Eq. ©, with T^{S) = 
4/3 « 1.33; S = {{a,d, f},{b,e,g},{c, f,h},{a,d,g},{b,e,h},{a,c, f},{b,d,g}, 
{c,e,h}) for Eq. ©, with T*^,{S) = 4/(8/3) = 1.5; and S = ({a, c, /}, e, g}, 
{c, e, /i}, {a, d, g}, {6, d, /, h}} for Eq. Q, with T*{S) = 4/(5/2) = 1.6. 

in S, there certainly always is a finite value of B, the number of buffer positions 
per node per path that goes through it, such that CI and C2 are satisfied. In 
all interleaved cases, this value is i? = 1. 

An example illustrating all of this is presented in Figure IH where we give a 
set of four source-to-destination paths, the graph G that eventually results from 
them, and also the three schedules that result in equalities in Eqs. ©-([S]). In 
this case the two inequalities in Eq. © are strict, since it can be shown that 
X(G) = 3, xr„t(G) = 8/3, and x*(G) = 5/2 [Ml [2]. 

5 Scheduling by edge reversal 

From the three schedules illustrated in Figure S] it would seem that finding 
a schedule S to maximize T{S) requires that we give up on the interleaved 
character of the underlying multicoloring and, along with it, give up on the 
equivalent property that edges of D that are consecutive on some source-to- 
destination path appear in S alternately. However, once color interleaving is 
assumed we are automatically provided with a principled way to heuristically try 
and maximize T{S) by appealing to a curious relationship that exists between 
multicolorings and the acyclic orientations of G. We now review this heuristic 
and later build on it by showing how to adapt it to abandon interleaving only 
on occasion during a schedule, aiming at maximizing T{S). 

An orientation of G is an assignment of directions to its edges. An orientation 
of G is acyclic if no directed cycles are formed. Every acyclic orientation has 
a set of sinks (nodes with no edges oriented outward), which by definition are 
not neighbors of one another. An acyclic orientation's set of sinks is then an 



10 



independent set. The heuristic we now describe, known as scheduhng by edge 
reversal (SER) [71 |5], is based on the fohowing property. Should an acyclic 
orientation be transformed into another by turning all its sinks into sources 
(nodes with no edges oriented inward) , and should this be repeated indefinitely, 
we would obtain an infinite sequence of independent sets, each given by the set of 
sinks of the current orientation. Though infinite, this sequence must necessarily 
reach a point from which a certain number of acyclic orientations gets repeated 
indefinitely. This follows from the facts that there are only finitely many acyclic 
orientations of G and that turning one of them into the next is a deterministic 
process. 

The orientations that participate in this cyclic repetition, henceforth called 
a period, have the property that every node of G appears as a sink in the same 
number of orientations. Furthermore, any two neighboring nodes of G are sinks 
in alternating orientations, regardless of whether the period has already been 
reached or not. It clearly follows that the sets of sinks in a period constitute 
a schedule that is based on an interleaved multicoloring. Depending on the 
very first acyclic orientation in the operation of SER more than one period 
can eventually be reached. The different periods' properties vary from one to 
another, but it can be shown that at least one of them corresponds to the optimal 
interleaved multicoloring, i.e., the one that yields xrnt(^) H- The heuristic 
nature of SER is then revealed by the need to determine an appropriate initial 
acyclic orientation. 

Determining a schedule S by SER follows the algorithm given next. We use 
ujQ, wi, . . . to denote the sequence of acyclic orientations of G. For t — 0,1, . . 
we denote by Sinks (wf) the set of sinks of ujt- 

Algorithm SER: 

1. Choose ojQ. 

2. t := 0. 

3. Obtain ujt+i from ujf 

4. If the period has not yet occurred, then t :— t + 1; go to Step 3. If it 
has, then let piojo) be its number of orientations, micjo) the number of 
times any node appears in them as a sink, and uJk,i^k+i, ■ ■ ■ j'^k+piujo)-! 
the orientations themselves. Output 

S = {S'mks{ujk), Sinks(a;fe+i), . . . , Sinks(wfe+p(„„)_i)) 

and 

T{S) = ^^^^^^ 

In this algorithm, the explicit dependency of both p(wo) and m(a;o) on ujq is 
meant to emphasize that, implicitly, the two quantities are already determined 
when in Step 1 the initial orientation cjq is chosen. As with the very existence of 
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the period, this follows from the fact that the algorithm's Step 3 is deterministic, 
so there really is no choice regarding the period to be reached once luq has been 
fixed. The role played by the two quantities is precisely to characterize the 
interleaved multicoloring mentioned above. That is, the period reached from 
Wo can be regarded as assigning q = m{uJo) distinct colors to each node of G 
using a total number p — p{i^o) of colors. Equivalently, it can be regarded as 
a schedule S for which delivered(iS) = Pq = Pm{u!o) (where the first equality 
is true of all multicoloring-based schedules, as we discussed in Section S]) and 
length(5) = L = p = p{ujo). By Eq. ([T]), the final determination of T{S) follows 
easily. 

As a final observation, we note that, although the knowledge of p(wo) and 
m{u!o) after Step 1 is only implicit, it can be shown that the ratio m(ajo)/p(wo) 
can be known explicitly at that point [7] . It might then seem that the remainder 
of the algorithm is useless, since the value of T{S) can be calculated right after 
Step 1. But the reason why the remaining steps are needed, of course, is that 

5 itself needs to be found, not just the T{S) that quantifies its performance. 

6 Improving on SER 

In Figure [5] we provide an example to illustrate why giving up interleaving may 
yield a schedule S of higher T{S). The general idea is that, given B, it may 
be possible to schedule a given transmission sooner than it normally would be 
scheduled by SER, provided there is a packet to be transmitted in the buffers of 
the sending node in D and the receiving node has an available buffer position 
for the path in question. While under SER two transmissions sharing a buffer 
alternate with each other in any schedule (and then B = 1 always suffices), 
disrupting this alternance implies that all buffering is to be managed in detail. 

In the example of Figure [5] transmissions g, h, and i are scheduled without 
regard to alternance if B = 2. While this results in improved performance (more 
packets delivered to node 4 per time slot), it is important to realize that this is in 
great part made possible by the structure of D in this example. Even though all 
three paths lead from node 1 to node 4, two of them are poised to interfere with 
each other particularly heavily by virtue of sharing node 2. The consequence 
of this is that transmissions on these paths will occur less in parallel than they 
might otherwise. But since B = 2 buffering positions are available per node 
per path, the path that goes through nodes 6 and 7 can compensate for this by 
transmitting twice as much traffic (thence the double occurrence of g in a row, 
and also of h and i, for each repetition of the schedule). This, however, is never 
detrimental to the traflac on the other two paths: all that is being done is to 
seize the opportunity to transmit in time slots that would otherwise go unused. 

Implementing the careful buffer management alluded to above requires that 
we look at the dynamics of acyclic-orientation transformation under SER in 
more detail. Given any acyclic orientation w of G, the node set of G can 
be partitioned into independent sets Ii, I2, ■ ■ ■ , Id such that Ii is the set of all 
sinks in G according to w, I2 is the set of all sinks we would obtain if all 
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Figure 5: A set of P = 3 directed paths (a), the resuhing directed multigraph D 
(b), and the resulting undirected graph G (c). The optimal SER schedule is <S = 
{{a],{h],{c,g},{d,i},{e,h},{f}), yielding T{S) = 3/6 = 0.5. An alternative 
schedule that does not comply with the SER alternance condition, with B = 2, is 
iS = {{g, c}, {g, /}, {h, &}, {h, e}, {i, a}, {i, d}), which results in an improvement 
to T{S) = 4/6 « 0.67. 

nodes in Ii were to be eliminated, and so on. In this partition, known as a 
sink decomposition, d is the number of nodes on a longest directed path of G 
according to lo. When oj is turned into lo' by SER a new sink decomposition is 
obtained, call it ![, I2, ■ ■ ■ , I'^, , such that ![ = I2, I2 =^ 6tc., with d' < d. The 
reason why equality need not hold in all cases, but set containment instead, is 
that each Ik may get enlarged by some of the previous orientation's sinks before 
becoming 

It is then possible to regard the operation of SER as simply a recipe for 
manipulating sink decompositions. At each iteration the set containing the 
current sinks is eliminated and its nodes are redistributed through the other sets. 
The remaining sets are renumbered by decrementing their subscripts by 1 and a 
new, greatest-subscript set may have to be created. The rule for redistributing 
each of the former sinks is to find the set of greatest subscript containing one 
of the node's neighbors in G, say 7^, and then place the node in set Ik+i- This 
is illustrated in panels (a) and (b) of Figure [51 

Altering this rule is the core of our modified SER, henceforth called SER 
with advancement (SERA). If i is the sink in question, the operation of SERA 
is based on placing node i in the least-subscript set that does not contain a 
neighbor of i in G. This clearly maintains acyclicity just as the previous rule 
does, but now the former sink is not necessarily turned into a source, but rather 
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Figure 6: Each set in a sink decomposition is represented by a rectangular box 
and numbered to indicate the set's subscript. Note that directed edges refer to 
acycHc orientations of G. Applying SER to the sink decomposition in panel (a) 
yields the one in panel (b). In this transformation both a and h are turned into 
sources. The alternative of using SERA, on the other hand, makes it possible 
for a to be placed in a lower-subscript set, avoiding the transformation into 
a source and yielding the sink decomposition in panel (c). This can be done 
only because the set to which a is added contains none of its neighbors in G. 
Assuming that transmissions a, e, and / are initially arranged consecutively in 
one of the P paths in the order e, a, f, we have i~ = e, i = a, and i+ = /. 
We also have, in reference to panel (a), = 4, = 2, and fc+ = 5. Then the 
additional conditions for the move from (a) to (c) to occur arc that the buffers 
shared by transmissions e and a contain at least one packet (since k < k~), and 
that those shared by a and / have room for at least one packet (since k < 



into a node that can have edges oriented both inward and outward by the current 

orientation, respectively from nodes in sets of greater subscripts and to nodes 
in sets of lesser subscripts. Additionally, while this alternative placement of 
node i does favor it by virtue of lowering the number of time slots that need 
to go by before it is a sink once again, clearly there is no detriment to any of 
the other transmissions, which will assuredly become sinks no later than they 
would otherwise. 

As we mentioned, however, unlike SER this rule can only be applied as 
a function of B and the buffering-related constraints we mentioned. Suppose 
that i is preceded by transmission i~ and succeeded by transmission i+ on the 
original path out of 1^1,1^2, ■■■ jVp to which it belongs. Of course, both i~ 
and i+ are nodes of G as well. Suppose further that these two nodes are in 
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sets /fe- and Jfc+, respectively, and that we are attempting to place node i in 
set Ik- The further constraints to be satisfied are the following. If fc < k~, 
then the buffers shared by transmissions i~ and i must contain at least one 
packet to be transmitted. If fc < then the buffers shared by transmissions 
i and i+ must contain room to store at least one more packet. This can all be 
implemented rather easily by keeping a dynamic record of all buffers. A simple 
case of evolving sink decompositions in the SERA style is shown in panels (a) 
and (c) of Figure [HI 

SERA, like SER, operates on finitely many possibilities and deterministically. 
A "possibility" is no longer simply an acyclic orientation, but instead an acyclic 
orientation together with a configuration of buffer occupation. In any event, 
periodic behavior is still guaranteed to occur and we go on denoting by p{uJo) 
the number of possibilities in the period that one reaches from ojq. The notion 
behind m(aJo), however, has been lost together with the certainty of interleaving, 
since SERA does not guarantee that every node of G is a sink in the period the 
same number of times. For i d N, an alternative definition is that of mi{uJo), 
which we henceforth use to denote the number of times node i is a sink in the 
period, not necessarily the same for all nodes. 

Determining the schedule S through SERA proceeds according to the fol- 
lowing algorithm. 



Algorithm SERA: 

1. Choose loq. 

2. t := 0. 

3. Obtain LUt+i from ujt, employing advancement as described. 

4. If the period has not yet occurred, then t := t + 1; go to Step 3. If it 
has, then let p{uJo) be its number of orientations (with associated buffer- 
occupation configurations), mi{u;o) the number of times node i appears in 
them as a sink, and uik,uJk+i, ■ ■ ■ ,u!k+p(u:o)-i the orientations themselves. 
Output 

S = (Sinks(a;fc), Sinks(a;fe+i), . . . , Sinks(wfc+p(^p)_i)) 

and 

T{S) = 



where T is the set of the nodes of G that correspond to terminal edges of 
the paths 'Pi, 7^2, • ■ • , Vp. 

In this algorithm, note that the determination of T{S) generalizes what is 
done in SER. This is achieved by adopting delivered(iS) = X^ieT '^i('^o) '^tiile 
maintaining length(5) = p(wo) in Eq. ([1]). Particularizing this to the case of 
SER yields delivered(iS) = Pm(cjo), as desired, since mi(wo) becomes m(wo) for 
any node i of G and moreover |r| = P. 
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7 Methods 



We have conducted extensive computational experiments to evaluate SER and 
SERA, the latter with a few different values for the buffering parameter B. Be- 
fore we present results in Section |S1 here we pause to introduce the methodology 
that was followed. This includes selecting the network topology that eventually 
leads to graph G and the choice of the initial acyclic orientation of G. 

7.1 Topology generation 

We generated 1600 networks by placing n nodes inside a square of side 1500. 
For each network the first node was positioned at the square's center. Given the 
nodes' communication (or interference) radius R, and with it the neighborhood 
relation among nodes (i.e., two nodes are neighbors of each other if and only 
if the Euclidean distance between them is no greater than R), we proceeded to 
positioning the remaining nodes randomly, one at a time. Positioning a node 
was subject to the constraints that it would have at least one neighbor, that 
no node would have more than A neighbors, and moreover that no two nodes 
would be closer to each other than 25 units of Euclidean distance. Repeated 
attempts at positioning nodes while satisfying these constraints were not allowed 
to number more than 1000 per network. When this limit was reached the 
growing network was wiped clean and a new one was started. The value of R 
was determined so that, had the nodes been positioned uniformly at random, 
a randomly chosen radius- i? circle would have expected density proportional 
to A/i?^ and about the same density as the whole network, i.e., A/i?^ oc n. 
Choosing the proportionality constant to yield R — 200 for n = 80 and A = 4 
results in the formula R = 200-y/20A/n. Of the 1600 networks thus generated, 
there are 100 networks for each combination of n G {60,80, 100, 120} and A G 
{4,8,16,32}. 

For each network we generated 50n sets of paths Pi,P2, ■ ■ ■ , Vp, each 100 
sets corresponding to a different value of P. Each of the sets resulted in a 
different D, then in G, as explained in Sections [5] and [3] The 50n sets comprise 
100 groups of n/2 sets each. The first of these sets for a group has P — 1 
and the single path it contains is the shortest path from a randomly chosen 
node to another in the network. Each new set in the group is the previous one 
enlarged by the addition of a new path, obtained by selecting two distinct nodes 
randomly, provided they do not already participate in the previous set. This 
goes on until P = n/2, so in the last set every one of the n nodes participates 
as either the origin or the destination of one of the P paths. For the sake of 
normalization, the results we present for T{S), given for P = 1, 2, . . . , n/2, are 
shown against the ratio P' = 2P/n G (0, 1]. 

7.2 Initial acyclic orientation 

Once G has been built for a fixed network and a fixed set of paths, the acyclic 
orientation ujq of G has to be determined. Our general approach is to label 
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every node of G with a different number and then to direct each edge from the 
node that has the higher number to the one that has the lower. Although the 
resulting orientation is clearly acyclic, we are left with the problem of labeling 
the nodes. We approach this problem by resorting to the paths Pi, P2, ■ • ■ , Vp 
from which G resulted, since the nodes of G are in one-to-one correspondence 
with the directed edges on the paths. It then suffices to number the paths' 
edges. 

We consider four numbering schemes: 

ND-BF. The paths are organized in the nondecreasing order of their numbers 
of edges (ties are broken by increasing path number). The edges are then 
numbered breadth-first from the path's origins, given this organization of 

the paths. 

ND-DF. The paths are organized in the nondecreasing order of their numbers 
of edges (ties are broken by increasing path number). The edges are then 
numbered depth-first from the paths' origins, given this organization of 
the paths. 

NI-BF. The paths are organized in the nonincreasing order of their numbers of 
edges (ties are broken by increasing path number). The edges are then 
numbered breadth-first from the paths' origins, given this organization of 
the paths. 

NI-DF. The paths are organized in the nonincreasing order of their numbers of 
edges (ties are broken by increasing path number). The edges are then 
numbered depth-first from the paths' origins, given this organization of 
the paths. 



8 Computational results 

We divide our results into two categories. First we give statistics on the 1600 
networks generated for evaluation of the algorithms. Then we report on the 
values obtained for T{S) by SER and SERA. 

One of the statistics is particularly useful: despite its simplicity, we have 
found it to correlate with the SERA results in a fairly direct way. This statistic 
is based on a function of G that aims to quantify how the interference among 
the initial P directed paths is reflected in the structure of G. We denote this 
function by p{G) and let it be such that 

In this equation, recall that the sets Yi, ^2, • • • , ^p, one for each of the initial 
directed paths, contain the edges that ultimately become the nodes of G. Thus, 
Y^p^i \^p\IP is the average number of edges on a path. Moreover, we let E' C E 
be the set of G's edges whose end nodes correspond to edges of distinct paths. In 
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words, then, p(G) is the average nuraber of off-path transmissions that interfere 
with the transmissions of a path having the average number of edges. 

8.1 Properties of the networks generated 

The 1600 networks' distributions of node degrees are given in Figure [71 which 
contains one panel for each of the four values of A and all four values of n. 
Their distributions of the numbers of edges on the P paths for P = n/2 are 
given in Figure |51 again with one panel for each of the four values of A and all 
four values of n. 

We see from Figure [7] that the degree distributions peak at the degree A, 
falling approximately linearly toward the lower degrees (except for A = 32, 
where a plateau is observed midway). Furthermore, the lowest observed degree 
grows with A, which is expected from the formula that gives the radius R as an 
increasing function of A. We also see from the figure that these distributions 
are approximately independent of the value of n for fixed A; in reference to that 
same formula, we see that letting R decrease with n does indeed have the ex- 
pected effect of maintaining an approximately uniform node density throughout 
the containing square. 

We also expect path sizes to be smaller as A increases, and this is in fact 
what Figure [8] shows. In fact, larger A values decreases the variability of path 
sizes, which moreover get concentrated around an ever smaller mean. For fixed 
A, what we see in the figure is a consistent shift to the right (i.e., greater path 
sizes) as n grows. This refiects the fact that larger n for fixed A leads to smaller 
i?, thus to longer paths. 

These observations are summarized in Table [U where the mean degree and 
mean path size are given for each combination of n and A values. This table 
also shows the average value of p{G), defined above as an indicator of how much 
interference there is in G among all P paths, when G refers to P = n/2. For 
fixed n, it is curious to observe that p{G) decreases as A is decreased from 32 
through 8, but then appears to fiatten out or even rebound slightly as A is 
further decreased to 4. Each of these averages corresponds to 10^ G instances 
(100 instances corresponding to the P = n/2 case of each of the 100 networks 
for fixed n and A) and is significant to the extent of the confidence interval 
reported for it in the table's rightmost column. As we demonstrate shortly, the 
peculiar behavior of p{G) helps explain a lot of what is observed with respect 
to how T{S) behaves in the case of SERA. 

8.2 Results 

Our results for SER are given in Figure |9] as plots of T{S) against the P' ratio 
introduced in Section [T] Each of the figure's four panels is specific to a fixed 
A value and shows a plot for each value of n combined with either the ND-BF 
or the ND-DF numbering scheme. All results relating to the NI-BF and NI- 
DF schemes are omitted, as we found them to be statistically indistinguishable 
from their ND counterparts. From this figure it seems clear that, as P' increases 
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Figure 7: Degree distributions of the 1600 networks, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). For each combination of n and A the distribution 
refers to 100 networks. 
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Figure 8: Path-size (number of edges) distributions of the 1600 networks, for 
A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). For each combination of 
n and A the distribution refers to 100 networks and to 100 path sets for each 
network, each set comprising P — n/2 paths. 
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Table 1: Mean values of the distributions in Figures [7] and HI and the average 
p{G) values for the 10^ G instances corresponding to each combination of n and 
A when P = n/2. Confidence intervals refer to these averages and are given at 
the 95% leveL 



n 


A 


Mean 


Mean 




P{G) 






degree 


path size 


Averag 


;e Conf. int. 


60 


4 


3.33 


7.46 


0.4 


0.06 




8 


6.22 


4.85 


0.37 


0.06 




16 


11.67 


3.57 


0.4 


0.05 




32 


21.23 


2.84 


0.58 


0.05 


80 


4 


3.36 


8.32 


0.37 


0.06 




8 


6.37 


5.36 


0.35 


0.06 




16 


12.17 


3.92 


0.39 


0.05 




32 


22.36 


3.06 


0.59 


0.05 


100 
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3.40 
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0.36 


0.06 




8 


6.40 


5.86 


0.34 


0.07 




16 
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4.22 


0.38 


0.06 




32 


23.09 


3.27 


0.6 


0.05 


120 


4 


3.40 


9.95 


0.34 


0.06 




8 


6.45 


6.28 


0.33 


0.07 




16 


12.50 


4.52 


0.38 


0.06 




32 


23.59 


3.47 


0.58 


0.05 



(i.e., as the number of paths P grows towards the superiority of the BF 

schemes over the DF schemes becomes apparent, more pronouncedly so for the 
lower values of A. The reason why the BF schemes tend to perform better 
than the DF schemes should be intuitively clear: the BF schemes number the 
transmissions that are closer to the paths' origins first, therefore with the lowest 
numbers. As the initial acyclic orientation of G is built from these numbers, 
the first sinks during the operation of SER will correspond to starting parallel 
traffic on as many paths as possible. Overall it also seems that larger values of 
n lead to better performance for fixed A, but the distinction appears to be only 
marginal and is sometimes obscured by the confidence intervals. 

A similar set of plots is given in Figure 1101 now displaying our results for 
SERA as plots of T{S) against the ratio P' . Once again there is one panel for 
each value of A, and once again several possibilities regarding the numbering 
schemes are omitted because of statistical indistinguishability. This is also true 
of the various possibilities for the value of B, with the single exception we 
mention shortly. Thus, most plots correspond to the ND-BF numbering scheme 
with B = 1. The single exception is that of n = 60 with A = 4, for which we 
also report on the B = 2 case. For fixed A and P', increasing n also leads to 
increased T{S). In the particular case of n — 60 and A = 4, increasing B from 
1 to 2 also causes T{S) to increase. 
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Figure 9: Behavior of T{S) for SER under the two numbering schemes ND-BF 
and ND-DF, with A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). Data 
are averages over the 10** G instances that correspond to each combination of n 
and A for each value of P. Error bars are based on confidence intervals at the 
95% level 
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Figure 10: Behavior of T{S) for SERA under the numbering scheme ND-BF, 
with A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over 
the 10^ G instances that correspond to each combination of n and A for each 
value of P. Error bars are based on confidence intervals at the 95% level. 
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As we fix A, n, and P', Figures and [TUl reveal tliat T{S) is consistently 
higher for SERA than it is for SER (by a factor of about 2 to 4) across all values 
for these quantities, thereby establishing the superiority of the former algorithm 
over the latter. For sufficiently large n this occurs for the same value of B (that 
is, for B = 1), which furthermore establishes that this superiority does not 
in general depend on the availability of more buffering space. It is, instead, 
determined solely by the elimination in SERA of the mandatory alternance of 
interfering transmissions in SER. 

Fixing n and P' while varying A (i.e., moving between panels) yields further 
interesting insight about the two algorithms. While for SER increasing A under 
these conditions causes T{S) to increase monotonically (though sometimes al- 
most imperceptibly) for the same numbering scheme, doing the same for SERA 
for constant B leads T{S) to behave in a markedly non-monotonic way. In fact, 
as A is increased from 4 to 8 there is also an increase in T{S), but increasing A 
further through A = 32 leads to decreases in T{S). As we anticipated earlier, 
this is fully analogous to the behavior of p{G) as A is increased in the same 
way while all else remains constant. This suggests that what determines the 
relative behavior of T{S) in these circumstances is the intensity of inter-path 
interference as it gets shaped by the structure of G. In other words, T{S) and 
p{G) tend to vary along somewhat inverse trends with respect to each other. 



9 Discussion 

When SER is used, it follows from our discussion in Sections S] and [5] that 
T{S) = T-l^iS). By Eq. ©, we then have 

where achieving equality requires that we choose optimally. Now let ip{G) = 
max{a;(G), \N\/a{G)}, where uj{G) is the number of nodes in the largest clique 
of G and a{G) is the number of nodes in the largest independent set of G. It 
can be shown that x*nt(G) > <yf(G)El 

whence 

P P'n 

^^'^ ^ Rg) = MG) ' ^''^ 

where we have taken into account the way we handle P in all our experiments. 
We see then that T{S) is bounded from above by the fraction of rt/2 given by 
P' /ip{G). For fixed P', this fraction tends to be small if the largest clique of 
G is large or its largest independent set is small, whichever is more influential 
on (f{G). Either possibility bespeaks the presence of considerable interference 
among the transmissions represented by the |A^| nodes of G. 



^See |18j . where the interleaved multichromatic number of G is referred to as G's circular 
chromatic number, and references therein. 
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Of course, in general we have no practical way of knowing how close each ujq 
we choose is to being the optimal starting point for SER, nor of knowing how 
different x*nti^) ^^'^ fi^) for the G instances we use. So the bound given in 
Eq. dTU]), located somewhere between 30P'/ip{G) and 60P'/(p{G) for our values 
of n, cannot be used as a guide to assessing how low the T{S) values shown in 
Figure IHl really are. But the bound's sensitivity to growing interference in G 
does provide some guidance, since all plots in the figure become flat from about 
P' — 0.3, regardless of the value of n or the numbering scheme used. Perhaps 
every G corresponding to such values of P' share some structural property, like a 
very large clique or only very small independent sets, that renders the resulting 
values of T(5) oblivious to all else. 

As for SERA, since the schedules it produces depart from a strict charac- 
terization as multicolorings of G, no upper bounds on T(S) are known to us. 
Nevertheless, a comparison with SER as provided by Figures M and [TO] reveals 
that T{S) for SERA surpasses T{S) for SER by a substantial margin, and also 
that SERA is capable of finding ways to improve T{S) somewhat even as P' 
grows. If our observation above regarding the structure of G as an inherent 
barrier to improving T{S) as P' grows is true, then the barrier's effects un- 
der SERA are considerably attenuated. This, we believe, is to be attributed 
to SERA's aggressively opportunistic approach of abandoning the interleaving 
that is the hallmark of SER. 

10 Concluding remarks 

Algorithms SER and SERA are methods for link scheduling in WMNs. As 
such, and unlike other methods for link scheduling, they are built around a 
set of origin-to-destination paths and aim to provide as much throughput on 
these paths as possible. From a mathematical perspective they are both related 
to providing the nodes of a graph with an efficient multicoloring, in the sense 
discussed in Section S) For SER this is strictly true, but for SERA the defining 
characteristic of a multicoloring, that each node receives the same number of 
colors, ceases to hold. As we demonstrated through our computational results 
in Section [HI it is precisely this deviation from the strict definition that allows 
SERA to surpass SER in terms of performance. 

The functioning of both SER and SERA is supported by the use of the inte- 
ger parameter B > 1, which indicates how many buffering positions each WMN 
node has to store in-transit packets for each of the paths that go through it. 
Choosing B = 1 suffices for SER because of its inherent property of alternating 
interfering transmissions, but B > 1 may in principle be needed for the advan- 
tages of SERA to become manifest. In the simulations we conducted, however, 
only rarely has this been the case, since on average increasing B beyond 1 pro- 
vided no distinguishable improvement. In this regard, we find it important for 
the reader to refer to Figure [5) once again. As we remarked upon discussing that 
figure, profiting from a, B > 1 situation under SERA is largely a matter of how 
uniformly interference gets distributed on the particular set of paths at hand. 



25 



Our results in Section[51 therefore, can safely be assumed to have stemmed from 
circumstances that, on average, led to highly uniformly distributed interference 
patterns. 

The centerpiece of both SER and SERA is the undirected graph G, which 
embodies a representation of all the interference affecting the various wireless 
links represented by the graph's nodes. As we explained in Section |31 the steps 
to building G depend on how one assumes the communication and interference 
radii to relate to each other, and also on which interference model is adopted. 
We have given results for a specific set of assumptions, but clearly there is 
nothing in either method precluding its use under any other assumptions: all 
that needs to be done is construct G accordingly. 

Analyzing either method mathematically is a difficult enterprise, but since 
their performance depends on the heuristic choice of an initial acyclic orientation 
of G, any effort profitably spent in that direction will be welcome. In addition 
to potentially better decisions regarding initial conditions, further mathematical 
knowledge on SER or SERA may also come to provide a deeper understanding 
of how upper bounds on T{S) relate to what is observed. As we mentioned 
in Section [9l one such bound is already known in the case of SER. Obtaining 
better bounds in this well as some bound in the case of SERA, remains 

open to further research. 

Another issue that is open to further investigation is how to handle the 
potential difficulties that SER may encounter in the face of a growing number 
of nodes in G [20l[T9]. These difficulties refer to the fact that, in the worst case, 
the time required to detect the occurrence of the period may grow exponentially 
with the square root of the number of nodes. They are inherited by SERA, since 
it generalizes SER, and may require the development of further heuristics if they 
pose a real problem in practice. In a related vein, sometimes it may be the case 
that only the value of T{S) is needed, not iS itself. Knowing the achievable 
throughput without requiring knowledge of the schedule itself can be useful for 
evaluating WMN topologies or routing algorithms for them. 

Should this be the case, then it is possible to estimate T{S) more efficiently 
than using the full-fledged algorithms we gave. We can do this in the case of 
SERA by recognizing that T{S) is the limit, as i ^> oo, of 

US) = ^'^-^7f°'^ (11) 

where TOi(tJo,i) is the total number of times node i appears as a sink in orien- 
tations ujo,ijJi, . . . ,ut. To see this, let o{t) denote any function of t such that 
limi^oo o{t)/t ^ 0. We then have X^^gr "^^('^o, = r{t) Y^ier mi{u}o) + o{t) and 
t + 1 = r{t)p{ujo) + o{t), where r{t) is the number of times the period has been 
repeated up to iteration t. The limit follows easily, and automatically holds 
also for SER by straightforward extension. The streamlined version of either 
algorithm consists simply of letting t evolve either through a sufficiently large 
value determined beforehand or until Tt{S) becomes stable. Any of the two 
alternatives does away with the need to detect the occurrence of the period. 
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We note, finally, that we have found the results given in Figure [TO] to be 
practically indistinguishable from those obtained through the strategy outlined 
above for the computation of T{S). We have verified this by letting T{S) = 
Tt+{S), where t+ is the least value of t for which \Tt{S) - Tt-US)\/Tt^wiS) < 
0.001. Here w is a window parameter and in our experiments we used w = \N\. 
As for this particular choice, it comes from realizing that in both SER and SERA 
it takes at most |iV| — 1 iterations for a node of G that is currently not a sink 
to become one. This, in turn, comes from the fact that in each iteration either 
algorithm necessarily decreases by 1 the number of edges on a longest directed 
path from any non-sink node to a sink. We can see that this is true of SER 
by viewing its dynamics in terms of how the orientations' sink decompositions 
evolve. We can see that it continues to hold in the case of SERA because SERA 
never places a former sink i into one of the sink-decomposition sets that already 
contains a neighbor of i in G (cf. Figure ^ . 
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